Vol. XX (XXXX) 



[PROC. XXXI SMP STYLE] 



No. X 



SHAPE-CHANGING COLLISIONS OF COUPLED BRIGHT SOLITONS 
IN BIREFRINGENT OPTICAL FIBERS 



M.Lakshmanan, T.Kanna and R.Radhakrishnan 
Centre for Nonlinear Dynamics 
Bharathidasan University 
Tiruchirapalli 620 024 
India 



We critically review the recent progress in understanding soliton propagation 
in birefringent optical fibers. By constructing the most general bright two-soliton 
solution of the integrable coupled nonlinear Schrodinger equation (Manakov model) 
we point out that solitons in birefringent fibers can in general change their shape 
after interaction due to a change in the intensity distribution among the modes 
eventhough the total energy is conserved. However the standard shape- preserving 
collision (elastic collision) property of the (1 + 1) dimensional solitons is recovered 
when restrictions are imposed on some of the soliton parameters. As a consequence 
the following further properties can be deduced using this shape-changing collision. 

(i) The exciting possibility of switching of solitons between orthogonally polarized 
modes of the birefringent fiber exists. 

(ii) When additional effects due to periodic rotation of birefringence axes are consid- 
ered, the shape changing collision can be used as a switch to suppress or to enhance 
the periodic intensity exchange between the orthogonally polarized modes. 

(iii) For ultra short optical soliton pulse propagation in non-Kerr media, from the 
governing equation an integrable system of coupled nonlinear Schrodinger equation 
with cubic-quintic terms is identified. It admits a nonlocal Poisson bracket structure. 

(iv) If we take the higher-order terms in the coupled nonlinear Schrodinger equation 
into account then their effect on the shape changing collision of solitons, during 
optical pulse propagation, can be studied by using a direct perturbational approach. 



1. Introduction 

The study of optical wave propagation in a nonlinear dispersive (dielectric) fiber has 
been receiving considerable attention in recent times as the fiber can support under 
suitable circumstances a stable pulse called optical soliton [1-3] . It arises essentially due 
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to a compensation of the effect of dispersion of the pulse by the nonhnear response of 
the medium. 

The analysis of such pulse propagation naturally starts from the Maxwell's equations 
for electromagnetic wave propagation in a dielectric medium, 



where the induced polarization for silica fibers is 



(1) 



(2) 



In equation(l), represents the electric field, c is the velocity of light, /xq and eo are the 
permeability and permitivity of free space respectively, x^™^ is the mth order suscepti- 
bility tensor. 

In order to analyse equation(l) it is necessary to make several simplifying assump- 
tions: (i) The nonlinear part of the induced polarization is treated as a small perturbation 
to the linear part, (ii) The optical field is assumed to maintain its polarization along the 
fiber, (iii) Fiber loss is assumed to be very small, (iv) The nonlinear response of the 
fiber is assumed to be instantaneous, (v) In a slowly varying envelope approximation for 
pulse propagation along the fiber, the electric field can be written as [4] 



F{x, y)E{z, t)e^^Koz-ujot) ^ ^ ^ 



(3) 



where c.c stands for complex conjugate, "e* is the unit polarization vector of the light 
assumed to be linearly polarized, E is the slowly varying electric field, F{x,y) is the 
mode distribution function in the {x, y) plane, while Kq and ojq denote the propagation 
constant and central frequency of the optical pulse respectively. 

Under the above assumptions, rewriting equation(l) by using the method of separa- 
tion of variables and introducing the coordinate system, T=t— — , moving with the pulse 

at the group- velocity Vg = one can obtain a wave equation for the evolution of E as 



dE K"d^E 



dz 



2 dT^ 



+ jo\E\'E = 0, 



(4) 



where 70= 



cA 



. Here A^f / denotes the effective core area of the single- mode fiber, n2 



1 



represents the nonlinear index coefficient and the parameter K" = 
accounts for the group velocity dispersion (GVD). After normalizing equation(4) and 



\K"\ 



using the transformation, q = 

z and r as t we get the ubiquitous nonlinear Schrodinger (NLS) equation, 

iq,-{sgnK")q„ + 2\qfq = 0, 



^ and then redefining ^ as 
) equation, 

(5) 
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in which Tq represents the width of the incident pulse, z and t are the normalized distance 
and time along the direction of propagation and q is the normalized envelope. The NLS 
equation admits the familiar bright envelope (one) soliton 



kiRCxp 


i {{kin 




cosh 




[t-kuz + v[°^) 





li^^t) = —r 7 ^^AT (6) 



in the anamolous dispersion region {K" < 0), where km and ku give the amplitude 
and velocity of the soliton respectively and rj^^ is an arbitrary complex parameter. In 
the normal dispersion regime {K" > 0), equation(5) possesses dark soliton solutions 
[5]. Equation(5) is valid in the picosecond regime since in obtaining the above equation 
the optical field is assumed to be quasi-monochromatic. For ultra short pulses (width 
< 100/s) one has to include additional effects such as third order dispersion, nonlinear 
dispersion, self induced Raman effect, etc. Then equation(5) is for example replaced by 
the higher-order NLS (HNLS) equation, 

iqz - {sgnK")qu + 2\q\^q - equt - Ai(|gpg)t - MqlU = 0. (7) 

Considerable attention has been paid to this HNLS equation which results from the 
delayed response of the fiber nonlinearity [6] . 

2. Electromagnetic Wave Propagation in Birefringent Fibers and the Cou- 
pled NLS Equation 

In general a single mode fiber can support two distinct modes of polarization which 
are orthogonal to each other. This phenomenon is known as birefringence. Among 
these two modes one corresponds to the ordinary ray (0-ray) in which the refractive 
index of the medium is constant along every direction of the incident ray. The other 
is the extraordinary ray (E-ray) whose refractive index for the medium varies with the 
direction of the incident ray. In an ideal fiber these two modes are degenerate, while in 
a real fiber due to the fiber nonlinearity this degeneracy is broken and the phenomenon 
is known as "modal birefringence" [1]. 

Thus due to the effect of birefringence and nonlinear response of the medium there 
is a possibility of interaction of two copropagating modes. As a result of this the phase 
of one mode not only depends on its own intensity (Self Phase Modulation (SPM)) but 
also on the intensity of the copropagating mode (Cross Phase Modulation (XPM)) [1,2]. 

The propagation equation for such modes can be obtained again from the Maxwell's 
equations by considering the electric field ^ in the slowly varying envelope approximation 
as 

^i^,t) = I [(e?i?i(z,t)e'^«^^ + e5i?2(z,t) e^^«^^)F(x, y)e-*'^''* + c.c] , (8) 

where Ei, E2 are the amplitudes of two polarization components and "e^i, "6*2 are the 
unit orthonormal polarization vectors and F{x, y) is the fiber mode distribution. 
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Proceeding in the same way as in the case of the single mode fiber one can obtain the 
following coupled system of equations for the envelopes of the two copropagating waves, 



I ^ K'' 
lEiz H Eit TT^ut 

Vgl 2 



iE2z H — 

Vg2 



K" 



E2u+lo{\E2? + B\Ei?)E2 = (i, 



(9) 



where z and t represent the normalized distance and time along the direction of propa- 
gation, Vgi and Vg2 represent the group velocity of the two copropagating waves Ei and 
E2 respectively and B = / ^® XPM coupling parameter (t? - birefringence 

ellipticity angle which varies between and J ). After suitable transformations as in the 
case of NLS equation we can rewrite the equation(9) in its normalized form as 



iqiz + gitt + 2M(|gi|2 + B\q2?)qi = 0, 
iq2z + q2tt + 2ii{\q2? + B\qi?)q2 =0. 



(10) 



The above equation is the coupled nonlinear Schrodinger (CNLS) equation which is 
in general nonintegrable. However for B=l this reduces to the celebrated Manakov 
equations [7], 



^qlz + q-^tt + '^^J■{Vh\ 



\q2?)qi 



0, 



ife + 92tt + 2/i(|5iP + |g2r)?2 = 0, 



(11) 



which is a completely integrable soliton system. The Lax pair for the Manakov system 
(11) can be identified as 



/ 




^{\<11? + \<12?)} 

-2Xql + iqtt 
-2Xq^ + iql^ 



{2Agi 


{2Ag2 \ 


+iqu} 


+iq2t} 


{2^2 


-iqtq2 






-iqiql 









(12) 



such that 



Lz-Mi + [L, M] = 0, 



(13) 



which is equivalent to the Manakov equations (11). The existence of infinite number of 
involutive integrals of motion confirms the integrability of the Manakov system (11). 



3. Bilinearization and Bright Two-Soliton Solution 

In recent years to study the solution properties of the integrable systems several 
effective tools have been developed which include inverse scattering transform method, 
Hirota's bilinearization method, Backlund transformation method, geometrical methods 
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and so on. In this section by applying Hirota's technique we point out that the most 
general bright one-soliton and two-soliton solutions for the Manakov system (11) can be 
obtained [8]. 

Considering equation (11) and by making the following bilinear transformation gi = 
f) 12 = J, where g(z,t), h(z,t) are complex functions while f(z,t) is a real function, the 
following bilinear equations can be obtained, 

{iD, + D^,)g.f = 0, {iD, + D^,)h.f = 0, D^f.f = 2^i{gg* + hh*), (14) 

where and Dt are Hirota's bilinear operators. The above set of equations can be 
solved by making the following power series exapansion to g, h and f: 

g = Xgi+X^g3 + ..., h = Xhi + X'ha + f = 1 + \^ f2 + X^f4 + (15) 

where A is the formal expansion parameter. The resulting equations, after collecting the 
terms with the same power in A, can be solved to obtain the forms of g, h and f. In order 
to get the one-soliton solution the power series expansions for g, h and f are terminated 
as follows, g=Xgi, h=Xh\ and / = 1 +A^/2. After following the procedure as mentioned 
before, the bright one-soliton solution is obtained as 



92 y \ P J 1 + e''i+''i*+« V ^2 y cosh(j7i7j + f)' 

where % = ^,-(^+^^+77^ , j = 1, = f , ^2 = 1, A = ^ + |/?|2), e« = 

^^(kl+k'Y ^ ' '^1 % ^'"^ complex parameters. Here ^/]l{Ai, A2) represents the 

unit polarization vector, km and k\i give the amplitude and velocity of the Manakov 
one soliton respectively. 

For obtaining the bright two-soliton solution, the series is terminated as g = Xg\ + 
A'^.gs, h = Xhi + A'^/is and / = 1 + A^/2 + A^/4. After proceeding in a similar fashion as 
in the case of the one-soliton solution, the following bright two-soliton solution with six 
arbitrary complex parameters fci, ^2, cti, a2 and /32 can be obtained, 

91 = ^ ^ 

^^gf/i ^ (32^^^ + e''^+''i'+''^+''i -|- e''^+''2+''2+'52 

92 = ^ , (17) 

where 

]j — I _^ ^m+vl+Ri _|_ gVi+vi+So _j_ gVl+m+Sg _|_ ^V2+vl+R2 _|_ Q-ni+ril+m+vl+Ri ^ 

50 ^ ^12 Ri ^ '^ll R2 ^ ^^22 

ki + k^' ki + fci* ' ^2 + ^2' 

51 ki-k2 I , 52 fc2 - fcl , . 

' = (fci + /fcr)(/fcr + /fc2)^"^"^^""^""^' ' =(fc2 + /fc2*)(fci + fc2)^"'""""^""^' 
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{ki + kt){kt + k2) 



(/3lK21 -/32Kll), 



^2 - fci 



|fcl-fc2p 



r(KiiK22 — K12K21) and K. 



{k2 + k^){ki + k^) 



(/32K12 - 



The above most general bright two-sohton solution corresponds to a shape changing 
(inelastic) collision of two bright solitons which will be explained in the following sections. 

4. Inelastic Collision and Switching of Bright Solitons in the Manakov Model 

The collision property of bright solitons can be revealed by analysing the asymptotic 
form of the two-soliton solution [8]. Without loss of generality we assume that kiR > 
and kii > fej, where i=l,2, which corresponds to a head-on collision. One can easily 
check that asymptotically the two-soliton solution becomes two well separated solitons 
Si and 52. Thus for the above condition, asymptotically the "qui's for the two-solitons 
become as (i)77i_R ~ 0, r]2R ±00 as z ^ ±00 and (ii)r?2fi ~ 0, tjir ^00 as z ^ ±00. 
This leads to the following asymptotic forms for the two-soliton solution. 
(i) Limit z ^ -00: (a) (ryii? ~ 0, r]2R -00) : 



91 

92 



ai 



-Ml- 



^2 



kiRc'^'^'sech riiR + 



Ri 



(18) 



where {Al~ ,Al^) — [fi{aial + Pi/3l)]~i {ai, Pi). Here superscript 1- denotes Si at the 
limit z ^ — 00 and subscripts 1 and 2 refer to the modes qi and q2. 
(b) 5*2 {r]2R ~ 77i_R ^ 00): 



12 



om+vi+Ra-Ri 



Al;- 



k2Re'''''sech {mR + 



{Rz - Ri) 
2 
(19) 

where (A?-,^^) = (^) c [^,ia2a*^ + 020^)]-^ [(ai,/3i)Kr/ - (a2,/32)K2i'] in which 



ai = {ki + k^)[{ki-k2){ala2+ (31132)]^ and c = 
(ii) Limit -z ^ 00: (a) Si {tjir ~0, r]2R 00): 

i 



Q2 



where {A\+,AI+) = (f|) c \pL{aial + pipl)]--^ [(ai,/3i)Kr2 - (a2,/32)K; 

a2 = {k2 + kt) [{ki - k2){aia*2 + pip^)]^ . 
(b) 52 {r]2R ~ 0, r]iR -00): 



92 



a2 
P2 



0V2 



2+ 



Ai 
Al+ 



where {AI+,AI+) = [M(a2a^ +/?2/32)]"'(a2,/32)• 



fc2i^e*''^^sech r]2R 



R2 
2 



{R3 - R2) 
2 

(20) 

n in which 



(21) 
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In the above set of equations, the solitons before interaction are given by equations 
(18) and (19) and the solitons after interaction arc given by equations (20) and (21). We 
observe from the above mentioned equations (18-21) that due to the interaction between 
two copropagating solitons Si and ^2, their amplitudes change from {Al^ km, A^^ km) 
and iAl~ k2R, Al~ k2R) to {Al^km,Al'^kiR) and {Al~^k2R, Af^k2R) respectively, in ad- 
dition to a change in their phases by an amount ^3-ih.-R2 g^,^^ -R3+R1+R2 respectively. 
The interesting behaviour which should be noted in this collision is that though the 
amplitude and phase of each soliton change during collision, the total intensity of each 
soliton is conserved, ie., -|- |A2''^p = i, where n=l,2 represent the solitons Si and 

5*2 respectively. Thus in such a collision there is a change in the distribution of intensity 
among the two component fields keeping the total intensity conserved. This is shown 
in Fig(l), where a head-on collision of two solitons is pictured for the parameter values, 
ki = 1 + i, k2 = 2 — 1, Pi = 1, ai = 1 and a2 = Here initially the time profiles 

of the two-solitons arc evenly split between the two components qi and (72. At the large 
positive z end the profile of the Si soliton is completely suppressed in the qi component 
while it is enhanced in the 92 component. Noticeable changes in S2 soliton also take 
place. 




Fig. 1: Intensity profiles and |g2|^ of the head-on collision of the two-soliton solution 
of the Manakov model. 

However we can regain the elastic collision for a2 = 1 [8]. Thus for the special case 
Q!i : Q!2 = /3i : /32 the standard elastic collision nature of the soliton can be obtained. The 
above analysis shows that the (H-l) dimensional soliton system given by equation(ll) 
exhibits a novel type of shape changing collision not seen in any other (1-1-1) dimensional 
evolution equation which led us to identify the exciting possibility of switching of solitons 
between modes by changing the phase [8] . This novel type of interaction led Jakubowski 
et al. to suggest a method for implementing computation in a bulk nonlinear medium 
without interconnecting discrete components [9]. 
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5. Periodically Twisted Birefringent Fibers and Soliton Interaction 

Soliton propagation in a periodically twisted birefringent fiber has gained considerable 
attention recently [1,10]. The effect of periodic twist on the soliton propagation can be 
described by the coupled wave equations of the form 

iqiz + qitt + Qqi+(^q2+'2.fi{\qi\'^ + \q2\'^)qi = 0, 

ife + g2« - ^92 + 0-gi + 2/x(|gi|^ + 192^)92 = 0, (22) 

where a and g are the normalized linear coupling constants caused by the periodic twist 
of the birefringence axes and the phase- velocity mismatch from resonance respectively. It 
can be easily verified that the transformation, qi = cosOe'^^^qiM — sm6e~^^^q2M, qi = 
s'mOe^^^qiM + cos6e~^^^q2M, where the subscript M refers to the Manakov model and 
r = _|_ (j2 g^jjjj — i^Q(^-i(£^^ reduces equation(22) to the coupled system given 

by equation(ll). Hence by making use of the solutions (16) and (17) we obtain the 
one-soliton and two-soliton solutions of equation(22) respectively as 

g'l = {cos^e'^'^^iM - sin^e~'^^A2M} ^c, 

92 = {sin^e'^MiM + cos^e-*^M2M}9c, (23) 

where qc = kiRexp{ir]ii)sech{riiR + ^) and 

[cos6'e*'^^Q!i - sin6'e-*r^/3i] e''^ + [cos6'e*i'^a2 - sm9e-'^^P2] e"^ 



91 



■ 



cos 



5)eir^+5i _sin6le-'r^+*i 



D 



+ ■ 



D 

cos 06^^''+^^ - sin6ie-'^^+*2l e''i+''2+')2 

, (24) 



D 

where all the parameters in (23) and (24) as well as the quantity D are defined in the 
previous section and q2 can be obtained by just replacing cos6 by sinO and sinO by 
— cos 9 in the above equation. The soliton interaction for this system can be studied by 
carrying out the asymptotic analysis as before. 

Then the form of and ^2 before interaction {z — oo) is given by 

(91,92) = {[cosOe^'Al- - sin^e-'^M^-], [smee^'Al' + cos^e-'^M^"]) q^' 

and 

(91,92) = ([cos6»e*rM2- - sin 6*6-*^^^^-], [sineie^^M^" + cosOe-'^'Al-]) q^' 

respectivcly, where q^~ = kijiexp{iriii)scc\i{riij^ + ^)., q^^ = k2Rexp{iri2i)scch{ri2R + 
{R3-R1) polarization vectors Aj^ , j=l,2 are the same as in equations(18,19). 
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Fig. 2: (a) Evolution of the intensity profiles and |g2p showing the suppression 
in switching between the two modes of soliton in a twisted birefringent fiber, (b) 
Evolution of the intensity profiles jg'ip and \q2\^ showing the suppression in switching of 
soliton and enhancement in ^2 soliton in a twisted birefringent fiber. 



Proceeding in a similar fashion the form of and ^2 after interaction (limit z— > oo) 

can be obtained as 

{qi,q2) = ([cos 6*6^^^ A J+ - smee-'^'Al+], [sin6le'^Mj+ + cos 9 e' Al+]) 

and 

(91,92) = {[cosee^^'A\+ - smee-'^'Al+], [smee^^' A\+ + cosee-^^' Al+]) 

respectively, where = kii^,exp{iriii)scc\i {r]iR + iii^i-^-^^ , = k2Rexp{ir]2i)sech{ri2ii+ 

^) and the polarization vectors ylj^,i,j=l,2 are defined in section-4. 

In order to facilitate the understanding of the above behaviour with reference to the 

optical soliton switching between the orthogonally polarized modes, it is convinient to 
obtain the oscillating part of the intensities associated with the above asymptotic forms as 

= \A\'^\^cosH+\A'J'^\^sin^0+i-iy\A'l'^\\A]^\sm20cos{2Tz + (^"^), l,j = 



= (2,t) 



1, 2, where (/ ^ j), ^"^ = tan~^ ( j - tan~^ ( — |^ j . Now let us analyse how the 



[Author and title] 



10 



presence of the oscillatory term in the above expression and the changes in A" affect 
the switching dynamics of the solitons and 5*2. 

case(i): For all the ^"^'s are non zero, there is a periodic intensity switching which 
is always present in both the solitons and in both the components before as well as after 
the interaction. 

case(ii): If any one of the ^"^'s is zero, and others are nonzero then soliton interac- 
tion suppresses or enhances switching dynamics. This is illustrated in Fig(2a), for the 
parameter values k-\_ = 1 + i, k2 = 2 — i, ai = (ii = p2 = 012 = p = 0.25 and 

a = 0.5. Here wc observe that for \A]^\ ^ the switching in the intensity of soliton Si 
is fully suppressed in both the modes gi and 52- 

casc(iii): If any two of the ^"^'s are zero then the soliton interaction suppresses 
and enhances switching dynamics. This can be verified from the Fig(2b), for the chosen 
parameters, ki = 1 + i(O.l), = 1 - i(O.l), ai = 0.86 + i(0.5), 02 = 0.5 + i(0.86), 
/Ji = 0.7 + i(0.72), /?2 = 0.44 + i(0.9) and p = cr = 0.25. Thus for this case due to the 
interaction there is an exchange of periodic intensity between the two modes of soliton 
S2 with the suppression of the switching dynamics in soliton Si. 

6. CNLS with Coupled Cubic-Quintic Nonlinearity 

The model for ultra short optical pulse propagation in non-Kerr media (in particular 
materials with high nonlinear coefficients even at moderate optical intensities, for exam- 
ple, semiconductor dopped glasses, organic polymers, etc.) can be obtained by expanding 
the induced polarization vector in Maxwell's equation(l) as 

'P{r,t)=eo[x^^l'^ + X^'^^'^^^ +X^'^^'^^^^^ + ■■■], (25) 

where x'"^^ is the fifth order nonlinear susceptibility and following the procedure as in the 
case of NLS equation. The resulting equation, in normalized form, describing the effects 
of quintic nonlinearity on the ultra short optical soliton pulse propagation is [11], 

iq, + qn + 2|g|\ + j\qfq + ijiqtti + n2{\qfq)t + h3{\qfq)t = 0. (26) 

The generalization of equation(26) in order to include midtimode effects leads to the 
following coupled cubic-quintic NLS equation[ll], after neglecting third order terms, 

iqi. + qitt + 2(|giP + |g2|')<Zi + (/OikiP + P2\q2\''fqi + 2p2[(Ti - Pi)ki|' 

+(t2 - P2)|92n |g2|^ qi - 2i[(pi|gi|2 + P2\q2\'^)qi]t + 2i{piqlqit + P2g2«2t)gi = 0, 

iq2z+q2tt + 2(|gi|^ + |g2|^)?2 + (nlgil^ +T2|g2|^)^g2 + 2ri[(/9i - ri)|gi|^ 

+{P2 - r2)\q2f] Wif 92 - 2i[{Ti\qif + T2\q2f)q2]t + 2i{Tiqlqit + T2q2q2t)q2 = 0, 

(27) 

where pi,p2,Ti and T2 are real free parameters. The integrable nature of the above 
equation can be studied by obtaining the Lax pairs and conserved quantities for it. The 
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Lax pair associated with equation(27) is 





M 





2\qi 


2\q2 \ 




+iqit 


+iq2t 


+ 192^) 


+eitqi 


+021 q2 




2iX^ 


-iqlq2 




-i\qi\' 






-iOiz 




-2\q2 


-iqiql 




+iq2t 




-i\q2? 



\ 



(28) 



^2*92 



-ie2z J 



where Oi = !^_^{pi\qi\'^ + P2\q2\'^)dt' , 6*2 = !*_^{Ti\qi\'^ +T2\q2?)dt' and A is the spectral 
parameter. The compatabiUty condition — Mt + [L,M] = leads to equation(27). 
Though the integrable system(27) possesses infinite number of conserved quantities, only 
the lower ones are of physical importance. They can be written as 



Cl = 



C2 = 



-- / dt[-i{qiql^ + q2q2t) + Pi\qi\^ + T2\q2\'^ 

* J —oo 



C3 



+{P2 + Ti)\qx\''\q2\% 

dt[{qiqlu + q2q*2u) + (kiP + 1921')' + i{\qi?Nu + \q2?N2t) 
(iV^kiP + 7V||g2|') 
12 and iVo = 



+2i{Niqxqlt + N2q2q*2t) 



(29) 



etc., where Ni = 6^ = Pikif + P2\q2\ and N2 = 021 = Tikil + T2m\- Here ci, C2 
and C3 may be related to the number operator, angular momentum and the Hamiltonian 
(energy) of the system(27) respectively. It is intriguing to note that the fields qa and 
g*, a=l,2 do not have canonical Poisson bracket relations. However under the following 
nonultralocal Poisson bracket structure, 

{qi{x),ql{y)} = 5{x-y) + ipie{x-y)qi{x)ql{y), 

{qi{x),qi{y)} = ipie{y - x)qi{x)qi{y), 

{qi{x),q2{y)} = -i{P20{x-y)-Tie{y-x))qi{x)q2{y), 

{qi{x),q^{y)} = i{p2e{x-y)-ne{y-x))qi{x)q*2{y), (30) 

etc., where €{x) = 9{x) — 0{—x) is the sign function defined through the step function, 
0{x) = 1 for x > 0, 9{x) = 0, for x < 0, the intcgrals(29) become involutivc. The 
interrelation between equation(27) and the Manakov system allows us to construct the 
soliton solution of equation(27) in terms of the known Manakov soliton solution. The fields 
of these two models are related through a nonlinear transformation for the dependent 
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iQlMt 




\Q2Mt\' 





(a) 



(b) 



Fig. 3: (a) Evolution of the intensity profiles jgiMiP and |g2MtP of the two-soliton 
solution of the Manakov model, (b) Intensity profiles \qit\'^ and |(72tP of the two-soliton 
solution of (27) for pi = p2 = n = T2 = 1. 



variables given by 



QaM = qaexp{-i9a), a = 1,2 



(31) 



and the subscrip M represents the Manakov model. Their Lax pairs are related through 
a local gauge transformation L' = g~^Lg — g~^gt, M' = g~^Mg — g~^gx where 



9 = 



10 

exp{-iei) 

exp{-ie2) 



(32) 



Then the one-soliton solution takes the form 

(91,92) = (ae'*i*^°^("(*-''^+*»,/Je'''^*^"^("(*-''^+^»)sec/i(iy(t - vz + 5))e'('**+'^^) (33) 

and the two-soliton solution can be obtained by substituting the two-soliton solution of 
the Manakov model in the transformation(31). 
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The effect of the dependence of phase change on the real free parameters pi, p2, n 
and T2 during colhsion can be revealed by carrying out the asymptotic analysis of the 
two-soliton solution as earlier. In this connection, we have taken the derivatives of the 
n,j=l,2, since the effect of the phases are reflected only in such terms rather than 

Wc have shown these effects for the parameter values fci = 1.5 + i(0.5), k2 = 2--i(0.7), 
cti = Pi = 132 = I and a.2 = '^'^gg^" , by comparing the plots of the intensity profiles in 
the absence of the free parameters (pi, p2, t\ and T2) with those in their presence. It 
is obvious from equation(31) that for the choice of parameters pi= P2=ti=T2= 0, the 
two-soliton solution of equation(27) reduces to that of the Manakov system. Hence in 
Figure(3a) we have plotted the intensity profiles of |giMtP and \q2Mt\^ for the above 
parametric values. In the figure there appears a splitting in each of the asymptotic 
profiles before and after interaction. 

After setting pi—p2—Ti=T2= 1, we can show that the splitting disappears as shown 
in Figure (3b). Thus in addition to the inelastic soliton collision, there is a change in the 
asymptotic form (suppression of splitting) of the intensity profiles which arises due to 
the presence of the phase terms Oa in the transformation(31). 

7. General Soliton Perturbation 

As mentioned in section(2) the CNLS equation is in general a nonintegrable one. We 
can apply the multiple scale perturbation theory to study such nonintegrable system with 
more general perturbations R\ and R2 as given below, 

i^iz + gitt + 2/x ( I gi I ^ + 1 92 1 ^ ) gi = i ei?i (gi , 92 , git , 92* • • • ) > 

ife + g2« + 2/i(|92|^ + |5ing'2 =ieR2{qi,q2,qit,q2t--)- (34) 

In the absence of the perturbation the above equation reduces to the Manakov system. 
After expanding qi as 

Qi = qf\kiR, hi, . *f ) + ^ + e'^r + - i, j=l,2, (35) 

where J^^^ and are the real and imaginary parts of the polarization vectors of the 
two-solitons respectively, in the asymptotic limits z ±00. Substitution of this in 
equation(34), it is trivial to verify that qf''' is nothing but the two-soliton solution of the 
Manakov system. Asymptotically this two-soliton solution of the Manakov model given 
by equation(17) can be written as a combination of two one-solitons as 2 — > ±00 

( q,\^( ql^+ql^\ 

By using a direct perturbational approach [12] the evolution equations for the soliton 
parameters in the presence of perturbations can be obtained. Applying this method to 
CNLS equation(lO) one can easily check that there occurs transmission and reflection 
scenarios during collision with a sensitive dependence of the collision outcome on the 



(36) 
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cross phase modulation coefficient and initial soliton parameters. As a special case by 
taking Ri = {—^Qi + (3\qj\'^qi + ^qut) where i,j=l,2 we have studied the effect of pertur- 
bations such as fiber loss, dispersion gain and nonlinear gain. The importance of taking 
such a form for perturbation is that in the absence of qi and 92 equation(34) along 
with the above specified form of Ri and R2 reduces to the Ginzburg-Landau equation 
governing pulse propagation in fiber amplifiers [1]. Since we have taken the asymptotic 
form of the more general two-soliton solution as the zeroth order solution, the collision 
properties of the system(34) can be directly studied by comparing the evolution of the 
soliton parameters in the limits z ^ —00 and z ^ 00. We expect that this may lead 
to some novel results such that collision among the solitons may be used to compensate 
fiber loss experienced by the interacting coupled one-solitons corresponding to any one 
of the polarized modes. 
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